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Abstract: Because the LHC is a proton-proton collider, sizable production of two pos- 
itively charged VF-bosons in association with two jets is possible. This process leads to 
a distinct signature of same sign high-p^ leptons, missing energy and jets. We compute 
the NLO QCD corrections to the QCD-mediated part of pp — > W + W + jj. These correc- 
tions reduce the dependence of the production cross-section on the renormalization and 
factorization scale to about ±10 percent. We find that a large number of W + W + jj events 
contain a relatively hard third jet. The presence of this jet should help to either pick up 
the W + W + jj signal or to reject it as an unwanted background. 
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1. Introduction 

The large energy and large luminosity of the LHC lead to sizable production cross-sections 
for multi-particle final states, including rather exotic ones. For example, one can produce 
a pair of positively charged 14^-bosons, in association with two jets pp — >• W + W + jj. At 
yfs = 14 TeV, the cross-section for this process is about 1 pb and therefore accessible. 
Leptonic decays of 14^-bosons give rise to two positively charged isolated leptons and missing 
energy, which is nearly a background-free signature. 

The observation of this process is interesting in its own right, but there are other 
reasons to study it. Of particular importance are various physics cases for which pp — > 
W + W + jj is a background process. Interestingly, such cases can be found both within 
and beyond the Standard Model. For example, it is possible to use same-sign lepton 
pairs to study double parton scattering at the LHC JIJ in which case the single scattering 
process pp —> W + W + jj is the background. Events with same-sign leptons, missing energy 
and two jets can also appear due to resonant slepton production which may occur in im- 
parity violating SUSY models J2] or in the case of diquark production || with subsequent 
decay of the diquark to e.g. pairs of top quarks. Similarly, one of the possible production 
mechanisms of the double-charged Higgs boson at the LHC has a signature of two same-sign 
leptons, missing energy and two jets Q|. 

At leading order in the perturbative expansion in QCD, the W + W + jj final state is pro- 
duced in proton collisions by both electroweak (EW) and QCD mechanisms. Interference 
terms between these mechanisms are doubly suppressed due to the different color structure 
of the two processes; interference only occurs at sub- leading color and even then only when 
the quarks are all identical. We neglect these terms and consider the QCD process sepa- 
rately. The corresponding cross-sections scale as <tew ~ a EW an< ^ ^QCD ~ «EW a s> wnere 
«ew and a s are electroweak and strong coupling constant respectively. Given the large 
hierarchy between strong and weak coupling constants, one expects ctqcd 3> cew> but this 



turns out to be too naive. In reality, the production cross-section due to a gluon exchange 
is only about fifty percent larger than the production cross-section of pp — > W + W + jj by 
electroweak mechanisms Q. 

As we mentioned earlier, W + W + jj production is the background to a number of 
interesting beyond the Standard Model (BSM) physics processes. It is peculiar that NLO 
QCD corrections to some of these signal BSM processes were calculated, see e.g. Refs. 0,13]. 
Also, the NLO QCD corrections to the electroweak production of W + W + jj were calculated 
recently in Ref. Q, but similar corrections to the QCD production of W + W + jj are 
unknown. There is a clear reason for this: the computation of NLO QCD corrections 
to the QCD-induced process pp — > W + W + jj is more involved because pp — > W + W + jj 
is a 2 — > 4 process. As a result, the NLO QCD calculation for QCD-induced pp — > 
W + W + jj requires dealing with one-loop six-point tensor integrals of relatively high rank, 
while for the EW-induced process this is not the case since at Born level there is no color 
exchange between the quark lines. Until very recently, theoretical methods for one-loop 
calculations were inadequate to handle computations of such a complexity, but the situation 
has changed dramatically in the past two years. Thanks to recent technical developments 
@, 0, |, |, |H| [11], 0, ||, 0, @, NLO QCD computations for such 2^4 processes as 



pp — > W(Z) + 3 jets, pp — )• tt + 66, qq — > 6666 and pp — > tt + 2 jets, have been performed 
during the past year pBl [L7], 18, 19, 20, 21, |2^, p3, p3, 25, 26|. Similar techniques should be 



applicable to NLO QCD computations of processes with several electroweak gauge bosons 
and jets, and pp — > W + W + jj is an interesting example. We also note that the computation 
of NLO QCD corrections to pp —> W + W + jj involves a small subset of amplitudes needed 
for the computation of NLO QCD corrections to pp — > W + W~jj, which is an important 
background to Higgs boson production in weak boson fusion. 

The remainder of the paper is organized as follows. In Section |2| we briefly discuss 
technical aspects of the calculation. In Section |3] we describe the results. In Section |] we 
present our conclusions. In the Appendix, numerical results for one-loop helicity amplitudes 
for — > {q~iW + qj)(q~kW + q m ) are given. 

2. Technical details 

In this Section, technical aspects of the computation are summarized. It is well-known 
that the computation of NLO QCD corrections to any process requires three ingredi- 
ents - one-loop virtual corrections to a Born process, real emission corrections and the 
subtraction counter-terms for infra-red and collinear singularities. To compute the one- 
loop virtual corrections, we use the framework of generalized Z)-dimensional unitarity 
|13, 14], closely following and extending the implementation described in Ref. [27]. Be- 



cause the one-loop virtual corrections are computed using unitarity cuts, an important 
ingredient for the computation of one-loop corrections are the tree-level helicity amplitudes 
— > (qiW + qj)(q~kW + q m ) + g. We compute those helicity amplitudes using Berends-Giele 
recursion relations p8| . Incidentally, these amplitudes are the ones that are needed to 
compute the real emission corrections to pp — > W + W + jj. The subtraction terms are cal- 



culated using the Catani-Seymour dipole formalism [29]. We employ the optimization of 
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Figure 1: A typical Feynman diagram that contributes to QCD production of W + W + jj in hadron 
collisions at leading order in perturbative QCD. 



the subtraction technique suggested in Ref. p0| ], where the subtraction is only performed 
if the kinematics of the event is close to either soft or collinear limit. Our implementation 



of the Catani-Seymour formalism closely follows the program MCFM [31|. In fact, we use 
the MCFM framework extensively to combine virtual and real emission corrections and the 
subtraction terms. Finally, we evaluate one-loop scalar integrals using the QCDIoop library 

H- 



We now describe the computation in some detail. First, we point out that within the 
generalized unitarity framework, a basic object that needs to be calculated for each phase- 
space point is a one-loop helicity amplitude. To decouple color, helicity amplitudes are 
expressed through primitive amplitudes [33|. Those primitive amplitudes can be computed 
with the help of the color-stripped Feynman rules |j^, ^, ^5|. We note that, for the 
calculation of a given primitive amplitude, color-charged particles are ordered while all 
permutations of color-neutral particles must be considered, to achieve a gauge-invariant 
result. For our purposes, this implies that the ordering of the W + bosons is not fixed and 
we have to account for all possible insertions of the W + bosons. Fortunately, given the fact 
that both of the two W bosons have the same charge, this is a relatively minor complication 
since two VF-bosons can not couple to the same quark line. Throughout the paper we treat 
the top quark as infinitely heavy, while all other quarks are taken to be massless. We do 
not consider final state top production as this leads to a different experimental signature. 
Top quarks therefore only contribute through fermion loops. While it is possible to deal 
with massive particles in the loop with D-dimensional unitarity [36, |37|, we choose to 
neglect these effects since the momentum transfer through the fermion loop is far below 
the threshold for top pair production. This is consistent with the treatment of parton 
distribution functions, which use /3q with five flavours in the evolution of a s . 



We now give the decomposition of the scattering amplitudes in terms of primitive 
amplitudes. We begin with the tree-level process — > udcs + {W + — > e + u e ) + {W + — > 
/x + f M ), where we fix the flavor structure, for definiteness and treat the Cabibbo-Kobayashi- 



Maskawa matrix as diagonal. The amplitude for this process reads 



A tTCC (u,d,c,s;e + ,u e ,fi + ,u IM ) = g 2 s I -y= j Pw+(s e+Ve )P w +(s^+ Vti ) 
I hdiAsiu ~ j-j-KiJisic) A (u,d,c,s), 



(2.1) 



where g s ,w are the strong and weak coupling constants, respectively, quark color indices 
are indicated explicitly, N c = 3 is the number of colors, and lepton labels in the right hand 



side of Eq^^) are suppressed. Also, we use the W + propagators with a Breit-Wigner 
form in Eq.( |2.1| ) 
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Finally, to account for leptonic decays W + — > Z + '(91)^7(92) i n Eq. fl2.lD , we replace the 
polarization vector of the outgoing W + bosons by 
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«(92)7m7-^(9i) 
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(2.3) 



(91+92) 2 ' ' 2 

A typical diagram that contributes to the amplitude Aq(u, d, c, s) is shown in Fig. [I]. 

The computation of real emission corrections requires the scattering amplitude for 
— > tides + g + (W + — > e + u e ) + (W + — > fi + v^). The amplitude for this process is written 
in terms of primitive amplitudes 

A tTee (ti,d,c,s,g;e + ,v e ,n + ,v^) = g z a I -j=j Pw+( s e+u e )Pw+( s ^+^) 
T isiu 5 idic A o(u, d, c, s, g) + T t a di J iai -A (ti, d, g, c, s) 



+ Jf { T i d iJisic A o(u,g,d,c,s)+T l a si _5 idi -A (u,d,c,g,s)) 



(2.4) 



Again, in the right hand side of this equation the lepton labels are suppressed. Finally, the 
one-loop amplitudes required for this calculation can be written as 

A°™- lo °v(u,d,c,s;e + ,v e ,v + ,v IM ) = 

9s ( ~j| ) p w+{s e + Ve )P w +{s^+ v ^) [5i d i- c 5 ia i^Ai + S id iJi a i 5 A 2 ] . (2.5) 

The color-ordered amplitudes A\^ can be expressed through the primitive amplitudes 
shown in Fig. 0. The relations read 

A\= ( N c - — ) A a (ti, d, c, s) - —A a (u, d, s, c) 



N r 



N r 



- —A b (ti, s, c, d) - —A c (u, s, c, d) + n f A d (ti, d, c, s) , 
A 2 = j-pA a (u, d, c, s) + I 1 + j-p J A a (u, d, s, c) 

+ Jp Ab ( U ' s ' 5 > d ) + 7^2 Ac ^ s > ^ d ) ~ jf Ad ( U > d ' 5 ' s "> 



(2.6) 
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Figure 2: Parent diagrams for one- loop primitive amplitudes for — > urfcs + W + W + . Shaded 
areas represent dummy lines which can not be cut. W + bosons are not shown. Quarks Q' do not 
couple to M^-bosons. 



Amplitudes with two adjacent quarks can be reduced to amplitudes where quarks and 
anti-quarks alternate using the C-parity relation 



A a (u Xu ,d Xd ,s Xs 



•) = -A a {u^,d Ad ,s A %c A °) 



(2.8) 



where we explicitly indicate fermion helicities 1 . 

Amplitudes with two identical quarks and anti-quarks are obtained from the amplitudes 
shown above by anti-symmetrising quarks or anti-quarks and including the appropriate 
symmetry factors in the cross-section calculation. Finally, we note that numerical results 
for helicity amplitudes at a particular phase-space point are given in the Appendix. 

3. Results 



In this Section, we present the results of the calculation. We consider proton-proton col- 
lisions with the center-of-mass energy y/s = 14 TeV. We require leptonic decays of the 
W-bosons and consider e + ^ + u e u^ final state. The VF-bosons are on the mass-shell and we 
neglect quark flavour mixing. We note that if identical leptons in the final state are present, 
there are interference effects not included in our calculation. However, such interference 
effects force the VK-bosons off the mass shell, so that their numerical importance is limited. 



1 We note that the r.h.s. of Eq. (2.S) is non-zero, since once the helicities of the fermions are fixed, we 
implement the coupling of the W-bosons to fermions as a vector coupling. 



Within this approximation, the cross-sections for same- and different-flavor production are 
related a(pp — > fj, + /j, + X) = a(pp — } e + e + X) = 0.5 a(pp — > e + fi + X). This implies that the 
cross-section for the full flavor sum e + fi + + e + e + + n + fi + can be obtained by multiplying 
our results by a factor two. 

We impose standard cuts on lepton transverse momenta p±j > 20 GeV, missing trans- 
verse momentum p_L im i ss > 30 GeV and charged lepton rapidity \r]i\ < 2.4. We define jets 
using anti-/cj_ algorithm [38|, with ARj ± j 2 = \/Jrjj[ — %>) 2 + {'Pji ~ 4>02) 2 = 0-4 and, unless 



otherwise specified, with a transverse momentum cut pj_j = 30 GeV on the two jets. The 
mass of the W-boson is taken to be mw = 80.419 GeV, the width 1~V = 2.140 GeV. W cou- 
plings to fermions are obtained from aQED(w-z) = 1/128.802 and sin 2 6\v = 0.2222. We use 
MSTW08LO parton distribution functions for leading order and MSTW08NLO for next- 
to-leading order computations, corresponding to a s (Mz) = 0.13939 and a s (Mz) = 0.12018 



respectively [39]. We do not impose lepton isolation cuts. All results discussed below apply 
to the QCD production pp — > W + W + jj; the electroweak contribution to this process is 
ignored. Note, however, that at next-to-leading order, QCD and electroweak production 
processes start to interfere and the separation of the two production mechanisms is not as 
clean as it is at leading order. 

We point out that, because a massive gauge boson is produced from each fermion line, 
the cross-section for the process pp — > W + W + jj remains finite even if the requirement 
that two jets are observed is lifted. Hence, we can consider the production of same-sign 
gauge bosons in association with n jets pp — > W + W + + n jets, where n = 0, 1, 2 or n > 2. 
It is useful to consider different jet multiplicities separately because, depending on the 
number of jets, backgrounds to pp — > W + W + +n jets change. Also, if we are interested in 
pp — > W + W + + n jets as a potential background to New Physics processes, it is helpful to 
know how jetty this process is. 

We begin by showing in Fig. |3| the dependence of the production cross-sections for 
pp — > e + / u + ^ e i/ M + n jets on the renormalization and factorization scales, which we set equal 
to each other. We show results for the following four processes i) pp — > W + W + + > 2 jets; 
ii) pp -> W + W + + 2 jets; Hi) pp -> W + W + + 1 jet; and iv) pp ->■ W + W + + jets. The 
total inclusive cross-section is then given by the sum of the two-jet inclusive cross-section, 
the one-jet exclusive cross-section and the zero-jet exclusive cross-section. The difference 
between inclusive and exclusive two jet cross-sections shows how often W + W + is produced 
in association with three, rather than two, jets. 

It is clear from Fig. g that leading order cross-sections for all jet multiplicities monoton- 
ically increase with the decrease of the renormalization/factorization scales. This behavior 
is driven by the dependence of the leading order cross-section on the square of the strong 
coupling constant; the factorization scale dependence is minor. It follows from Fig. || 
that when cross-sections are calculated at next-to- leading order in perturbative QCD, they 
show significantly reduced scale dependence. For example, considering the range of scales 
50 GeV < fj, < 400 GeV, we find the two-jet inclusive cross-section to be <t lo = 2.7 ± 1.0 fb 
at leading order and <t nlo = 2.44 ±0.18 fb at next-to-leading order. The forty percent 
scale uncertainty at leading order is reduced to less than ten percent at NLO. We observe 
similar stabilization of the scale dependence for the 0- and 1-jet exclusive multiplicities. 
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Figure 3: The dependence on factorization and renormalization scales of cross-sections for pp — > 
e + fi + v e Vp, + n jets, n = 0, 1, 2 at leading and next-to-leading order in perturbative QCD We set 
the two scales equal to each other /ip = /iR = /i. 



Combining these cross-sections we obtain a total NLO cross-section of about 2.90 fb for 
pp — > e + /i + z/ e i/ M inclusive production. This implies about 60 e + /j, + + e + e + + fi + fj> + events 
per year at the LHC with 10 fb~ annual luminosity. While this is not a gigantic number, 
such events will have a very distinct signature, so they will definitely be seen and it will 
be possible to study them. 

On the other hand, it is apparent from Fig. Blthat there is quite a dramatic change in the 
two-jet exclusive cross-section. At leading order, there is no difference between the exclusive 
and inclusive cross-sections, but this difference appears at NLO. With ARjj = 0.4 and a 
transverse momentum jet cut of 30 GeV, the two-jet exclusive cross-section is only about 
0.7 fb at ii = 150 GeV; if the scale is decreased to about 80 GeV, the two-jet exclusive cross- 
section becomes negative. One can argue that this is the consequence of the fact that the 
30 GeV jet transverse momentum cut is too small for the convergence of the perturbative 
expansion of the two-jet exclusive cross-section. However, it is not fully clear how to make 
this explanation compatible with the reasonable perturbative stability of the zero-jet and 
the one-jet exclusive cross-sections for the same value of the jet transverse momentum cut. 
In Fig. 0] we show the inclusive and exclusive cross-sections for pp — > e + fi + u e i/ fl + 2 jets at 
next-to-leading order in perturbative QCD, in dependence on the jet p± cut. We set the 
factorization and renormalization scales equal to each other [ir = fip = /u, and vary the 
scale [i in the range 100 GeV < /i < 200 GeV. It is clear that for a jet cut of 40 — 50 GeV, 
the scale dependence of the exclusive two-jet cross-section is relatively small, while for 
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Figure 4: The dependence on the jet p± cut of the two-jet inclusive and two-jet exclusive cross- 
sections for pp — > e + p + v e V/j, + 2 jets at leading and next-to-leading order in perturbative QCD. 
We show scale uncertainty bands, for values of the renormalization and factorization scales set to a 
common value /i, which is varied in the interval 100 GeV < p < 200 GeV. Results for fi = 140 GeV 
arc shown as solid lines. 



the 30 GeV jet cut, the scale dependence increases strongly. The results displayed in 
Fig. f| suggest that, whatever the exact value of the exclusive two-jet cross-section is, it is 
significantly smaller than the two-jet inclusive cross-section. This smallness implies that 
quite a large fraction of events in pp — > e + fi + v e v lJL + > 2 jets have a relatively hard third jet. 
This effect remains present, although less pronounced, for larger values of R. This feature 
may be useful for rejecting contributions of pp — > W + W + jj when looking for multiple 
parton scattering. 

We now turn to the discussion of kinematic distributions. Unless explicitly stated 
otherwise, we will consider two-jet inclusive processes. A glance at Fig. || suggests that 
NLO QCD corrections to jet cross-sections are small if the renormalization and factorization 
scales are set to /i = 150 GeV. Indeed, for /x = 150 GeV, two-jet inclusive leading and 
next-to-leading order cross-sections nearly coincide ctlo = 2.4 fb and onlo = 2.5 fb, so that 
the NLO QCD corrections change the leading order cross-sections by less than five percent. 
In the plots that follow, we show the scale uncertainty bands for 50 GeV < \x < 400 GeV 
and predictions for /i = 150 GeV, for leading and next-to-leading order distributions. 

The transverse momentum distributions of the hardest and next-to-hardest jets, as 
well as the distribution of the scalar sum of the transverse momenta of all jets and charged 
leptons and the missing transverse momentum in the event are shown in the first three pane 
of Fig. ^. It is clear from Fig. [| that jets in pp — > W + W + jj are hard; a typical transverse 
momentum of the hardest jet is close to 100 GeV and the transverse momentum of the 
next-to-hardest jet is close to 40 GeV. While no major changes in kinematic distributions 
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Figure 5: Kinematic distributions of the two hardest jets in the process pp —¥ e + /i + v e v M + 2 jets 
at leading and next-to-leading order in perturbative QCD for inclusive two-jet events. The bands 
show renormalization and factorization scale uncertainty, for 50 GeV < /x < 400 GeV. Solid 
lines show leading and next-to-leading order predictions for \x = 150 GeV. We use -Ht,tot = 
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occur when NLO QCD corrections are calculated, some shape changes are apparent from 
Fig. ||. Indeed, the NLO QCD transverse momenta distributions of the hardest and next- 
to-hardest jet show a characteristic depletion at large values of p±j- One reason this 
change occurs is because a constant, rather than a dynamical, renormalization scale is 
used in our leading order calculation. As was emphasized many times recently, the choice 



of dynamical re-normalization scales in leading order computations may affect shapes of 



kinematic distributions in a manner similar to NLO QCD corrections ||, 2C, 23, [Kj. The 
shapes of rapidity distributions of the hardest and next-to-hardest jets, as well as the 
distribution of the rapidity differences between hardest and next-to-hardest jets, shown in 
the next three pane, remain very similar to the shapes of leading order distributions, once 
the NLO QCD corrections are included. On the other hand, scale dependencies of rapidity 
distributions are reduced dramatically. 

We now turn to lepton and missing transverse momentum distributions. In Fig. pi 
we show the charged lepton transverse momentum distribution, the missing transverse 
momentum distribution, the distributions of the charged lepton rapidity and of the rapidity 
difference of the two leptons, as well as the distributions of the invariant mass of the 
two charged leptons and the transverse mass of the two VF-bosons 2 . Lepton kinematic 
distributions are affected by the QCD radiative corrections in a similar way to the jet 
kinematic distributions. We observe that transverse momentum distributions become softer 
while shapes of rapidity distributions are unaffected. Distributions of dilepton invariant 
mass and the transverse mass of the VF-bosons become softer as well. It is also interesting 
to look at the angular distributions of the charged leptons. In Fig. [?] the angular distance 
ARij = ^J{i]i — rjj) 2 + {<pi — (fij) 2 between a charged lepton of fixed flavor (e + or /x + ) and 
the hardest (next-to-hardest) jet is displayed, as well as the distribution of the relative 
azimuthal angle of the two charged leptons. Although distributions of angular distances 
between leptons and jets are broad, they peak at ARij ss 3 for both hardest and next- 
to-hardest jets. NLO QCD effects do not change this conclusion but, interestingly, they 
make the angular distance between next-to-hardest jet and the charged lepton somewhat 
larger at next-to- leading order. The distribution of the relative azimuthal angle of the two 
charged leptons becomes less peaked at A<^+j+ = n, although the two leptons still prefer to 
be back to back. It is interesting to remark that, if the two same sign leptons are produced 
through a double-parton scattering mechanism, their directions are not correlated. Hence, 
yet another possibility to reduce the single-scattering-background is to cut on the relative 
azimuthal angle between the two leptons. 

4. Conclusion 

In this paper, we presented the computation of NLO QCD corrections to the QCD-mediated 
process pp — > W + W + jj at the LHC. In spite of the fact that this is a 2 — > 4 process, the 
computation of NLO QCD corrections is relatively straightforward, thanks to spectacular 
developments in technology of NLO QCD computations in recent years || [7|, ||, ^, |l0, 11, 12 



13, 14, |lq|. We note that NLO QCD corrections to electroweak-mediated pp — > W + W + jj 



process were recently calculated in Ref. ||. 

2 Since neutrinos are not observed, it is impossible to reconstruct the true transverse mass of the pair of 
the two W-bosons. We follow Ref.Q and define m 2 ±ww — (E ±tl+l+ + -Ex.miss) 2 — (p±,i+i+ + Px,miss) 2 , 
where the missing transverse energy -Ex, miss is reconstructed from the missing transverse momentum using 
the invariant mass of the charged lepton system -Ex, miss = \/Px m i ss + m f+i- ■ 
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Figure 6: Kinematic distributions of charged leptons and missing energy in pp — > e + /i + v e v^ + 
n jets at leading and next-to-leading order in perturbative QCD for inclusive two-jet events. The 
bands show renormalization and factorization scale uncertainty, for 50 GeV < /i < 400 GeV. Solid 
lines show leading and next-to-leading order predictions for /i = 150 GeV. 



The production of two same-sign W-bosons in association with two jets is a background 
to studies of double-parton scattering, as well as to a number of beyond the Standard Model 
physics processes. The NLO QCD corrections to pp — > W + W + jj reduce the scale depen- 
dence to about ten percent. An interesting feature of pp — > W + W + jj is that perturbative 
QCD is applicable even if the observation of the two jets is not demanded. This opens up 
the possibility to study QCD radiative corrections to pp — > W + W + +n jets where the num- 
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Figure 7: Angular kinematic distributions in the process pp —> e + /i + ^ e ^ M + 2 jets at leading and 
next-to-leading order in perturbativc QCD for inclusive two-jet events. The bands show renormal- 
ization and factorization scale uncertainty, for 50 GeV < /i < 400 GeV. Solid lines show leading 
and next-to-leading order predictions for /i = 150 GeV. 



ber of jets can be zero, one and two. At next-to-leading order, the final state with three jets 
also appears. We have studied QCD corrections to processes with various jet multiplicities 
and observed reasonable behavior of the perturbative expansion for two-jet inclusive, and 
one- and zero-jet exclusive processes. The situation with the two-jet exclusive cross-section 
is less satisfactory. As follows from Fig. ||, the two-jet exclusive cross-section is about fifty 
percent of the inclusive cross-section. Since at leading order exclusive and inclusive two-jet 
cross-sections coincide, the NLO QCD correction is large. In addition, there is significant 
residual scale uncertainty in the prediction for the exclusive two-jet cross-section which is 
smaller than, but comparable to, the leading order scale dependence. We believe that these 
results suggest that nearly fifty percent of all events in pp — > e + fi + u e u fM + > 2 jets contain 
a relatively hard third jet. It may be possible to use this effect in designing selection or 
suppression criteria for W + W + jj final state. We have also studied a variety of kinematic 
distributions and showed that no dramatic shape changes occur. 

As a final comment, we point out that the production of two negatively charged W- 
bosons at the LHC is also possible, at a rate which is about forty percent of the W + W + 
production O. We can use our set up to compute NLO QCD corrections to pp — > W~W~jj 
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using charge conjugation and parity reversal ||. To this end, if we are interested in 
pp — > e~fi~i> e r> fl jj, we can do a calculation for pp — > e + fi + u e i , fl jj, treat the final state 
leptons as if they are positively charged and reverse the momentum directions for parity- 
odd distributions. 
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A. Results for primitive amplitudes at a fixed phase space point 

In this Appendix we give numerical results for all primitive amplitudes needed to recon- 
struct the one-loop amplitudes in Eq. (2.5) at a particular phase space point. The momenta 



(E,p x ,Py,p z ), in GeV, are chosen to be 

k^u) = (-500.000000000000, -500.000000000000,0.00000000000000,0.00000000000000), 
k 2 (d) = (-500.000000000000, 500.000000000000, 0.00000000000000, 0.00000000000000), 
fcg(c) = (54.2314070117999, -31.1330162081798, -7.92796656791140, 43.6912823611163), 
k 4 ( s ) = (214.488870161418, -27.0607980217775, -98.5198083786150, 188.592247959949). 
k 5 (e + ) = (85.5312248384887, -8.22193223977868, 36.1637837682033, -77.0725048002414), 
k 6 (v e ) = (181.428811610043, -57.8599829481937, -171.863734086635, -5.61185898481311), 
k 7 (n + ) = (82.8493010774356, -65.9095476235891, -49.8952157196287, 5.51413360058664), 
fc 8 (i/ M ) = (381.470385300815, 190.185277041519,292.042940984587, -155.113300136598). 

(A.l) 

Because the VF-bosons couple only to left-handed fermions and right-handed anti- 
fermions, all helicities in the following are fixed, so that {u, d, c, s, e + , z^ e ,/i + , u^\ have 
helicities {+,—,+, —,+,—,+, —}. Below we tabulate numerical results for the primitive 
tree and one-loop amplitudes at the phase-space point given in Eq. (A.l). In Table || we 



display tree-level amplitudes and ratios of unrenormalized one-loop amplitudes to tree- 
amplitudes 

n = ~^ (A.2) 



where cr = ^.^ /t —^ and the Ai are defined through equations ( |2.5[ ) - (^7J). 



r(i + e)r 2 (i-e) 

(4vr) 2 - e r(l - 2e) 
We use hr = 80 GeV; it enters our calculation through the usual modification of the 

loop integration measure J d<* -> „| J i"k in dimensional regulation. The one-loop 

amplitudes are calculated in the four-dimensional helicity scheme [42]. In Table y we give 

the ratio S 

4?r Re(A one - lo °P^ tree *) 

a s \A tree2 \ y ' 

of amplitudes summed over spin and color. Note that all unrenormalized virtual amplitudes 
(including fermion loops) contribute to yl one - lo °P. 



14 



Amplitude 


iA 2 


1/e 


e° 


Ao(u,d,c,s) 






-7.488599 - i 17.47659 


r a (u,d,c,s) 


-2.000000 


5.580909 + i 0.000000 


11.010146 + i 14.39329 


r a (u,d,s,c) 


-2.000000 


6.634395 + i 0.000000 


-16.15227-^7.909865 


r b (u,s,c,d) 


-1.000000 


3.551457-^3.141593 


-2.605243 + i 9. 184110 


r c (u,s,c,d) 


-1.000000 


-2.092934-i3.141593 


-1.237780 -i 9.355851 


r d (u,s,c,d) 




-0.666667 + i 0.000000 


1.431296 -i 2.305401 



Table 1: Numerical results for the primitive tree- level amplitude Ao(u, d,c,s), in units of 
1CP 11 GcV~ and the ratios of primitive one-loop amplitudes r;. 



Ratio 


1/e* 


Ve 


e° 


S(u,d, c, s) 


-5.333333 


13.62554 


23.35965 



Table 2: Numerical results for the ratio of spin and color summed amplitudes S(u,d,c,s). 

References 

[1] A. Kulcsza and W. J. Stirling, Phys. Lett. B475, 168 (2000); 
E. Maina, JHEP 0909, 081 (2009); 
J. R. Gaunt, C. H. Kom, A. Kulesza and W. J. Stirling, hep-ph/1003.3953. 



[9 
[10 

[11 

[12 
[13 

[14 
[15 

[16 



H. K. Dreiner, S. Grab, M. Kramer and M. K. Trenkel, Phys. Rev. D75, 035003 (2007). 

T. Han, I. Lewis and T. McElmurry, JHEP 1001, 123 (2010). 

See e.g. J. Maalampi and N. Romanenko, Phys. Lett. B 532, 202 (2002) and references 
therein. 

B. Jager, C. Oleari and D. Zeppcnfeld, Phys. Rev. D80, 034022 (2009). 
A. Dcnncr and S. Dittmaier, Nucl. Phys. B 734, 62 (2006). 

T. Binoth, J. P. Guillct, G. Heinrich, E. Pilon and C. Schubert, JHEP 0510, 015 (2005). 

R. Britto, F. Cachazo and B. Feng, Nucl. Phys. B 725, 275 (2005). 

R. Britto, F. Cachazo and B. Feng, Phys. Lett. B 611, 167 (2005). 

D. Forde, Phys. Rev. D 75, 125019 (2007). 

G. Ossola, C. G. Papadopoulos and R. Pittau, Nucl. Phys. B 763, 147 (2007). 

G. Ossola, C. G. Papadopoulos and R. Pittau, JHEP 0805, 004 (2008). 

R. K. Ellis, W. T. Giele and Z. Kunszt, JHEP 0803, 003 (2008). 

W. T. Giele, Z. Kunszt and K. Mclnikov, JHEP 0804, 049 (2008). 

C. F. Bergcr et al., Phys. Rev. D 78, 036003 (2008). 

A. Bredenstein, A. Denner, S. Dittmaier and S. Pozzorini, Phys. Rev. Lett. 103, 012002 
(2009). 



[17] A. Bredenstein, A. Denner, S. Dittmaier and S. Pozzorini, JHEP 1003, 021 (2010). 



15 



[18] G. Bevilacqua, M. Czakon, C. G. Papadopoulos, R. Pittau and M. Worck, JHEP 0909, 109 
(2009). 

[19] C. F. Berger et al., Phys. Rev. Lett. 102, 222001 (2009). 

[20] C. F. Berger et al., Phys. Rev. D 80, 074036 (2009). 

[21] R. K. Ellis, K. Mclnikov and G. Zanderighi, JHEP 0904, 077 (2009). 

[22] R. Keith Ellis, K. Melnikov and G. Zanderighi, Phys. Rev. D 80, 094002 (2009). 

[23] K. Mclnikov and G. Zanderighi, Phys. Rev. D 81, 074025 (2010). 

[24] T. Binoth, N. Greiner, A. Guffanti, J.-Ph. Guillet, T. Reiter, J. Rcutcr, Phys. Lett. B 685, 
293, (2010). 

[25] G. Bevilacqua, M. Czakon, C. G. Papadopoulos and M. Worek, Phys. Rev. Lett. 104, 162002 
(2010). 

[26] C. F. Berger et al., arXiv: 1004. 1659 [hep-ph]. 

[27] R. K. Ellis, W. T. Giele, Z. Kunszt, K. Melnikov and G. Zanderighi, JHEP 0901, 012 (2009). 

[28] F.A. Bcrcnds and W. T. Giele, Nucl. Phys. B 306, 759 (1988). 

[29] S. Catani and M. H. Seymour, Nucl. Phys. B485, 291 (1997) [Erratum-ibid. B 510, 503; 
(1998)]. 

[30] Z. Nagy and Z. Trocsanyi, Phys. Rev. D 59, 014020 (1999) [Erratum-ibid. D 62, 099902 
(2000)]; Z. Nagy, Phys. Rev. D 68, 094002 (2003). 

[31] J. M. Campbell and R.K. Ellis, Phys. Rev. D 62, 114012 (2000). The MCFM program is 
publicly available from http://mcfm.fnal.gov. 

[32] R.K. Ellis and G. Zanderighi, JHEP 0802, 002 (2008). 

[33] Z. Bern, L. J. Dixon and D. A. Kosowcr, Nucl. Phys. B 437, 259 (1995). 

[34] Z. Bern, L. J. Dixon and D. A. Kosowcr, Nucl. Phys. B 513, 3 (1998). 

[35] Z. Bern, L. Dixon and D. Kosower, Ann. Rev. Nucl. Part. Sci. 46, 109 (1996). 

[36] R. K. Ellis, W. .T Giele, Z. Kunszt and K. Melnikov, Nucl. Phys. B 822, 270 (2009). 

[37] K. Mclnikov and M. Schulze, Nucl. Phys. B 840, 129 (2010). 

[38] M. Cacciari, G. P. Salam and G. Soycz, JHEP 0804, 063 (2008). 

[39] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63, 189 (2009). 

[40] Contribution by S. Hoche, J. Huston, D. Maitre, J. Winter and G. Zanderighi to T. Binoth et 
al. "The SM and NLO multileg working group: Summary report", arXiv: 1003. 1241 [hep-ph]. 

[41] A. Alves, O. Eboli, T. Plehn and D. Rainwater. Phys. Rev. D 69, 075005 (2004). 

[42] Z. Bern and D. A. Kosower, Nucl. Phys. B 379, 451 (1992); Z. Bern, A. De Freitas, 
L. J. Dixon and H. L. Wong, Phys. Rev. D 66, 085002 (2002). 



16 



